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Abstract 

The CompHEP software package allows the evaluation of cross section and decays of elementary 
particles with a high level of automation. Arbitrary tree level processes can be calculated starting 
from the set of vertices prescribed by a given physical model. This note describes the details of the 
Minimal Supersymmetric Standard Model (MSSM) implementation in the CompHEP package, and 
the notation for the particles and parameters of the MSSM in CompHEP. 
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1 Introduction 
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The CompHEP package Q allows to generate automatically the Feynman diagrams for a process describ- 
ing collisions or decays of elementary particles starting from the Feynman rules. CompHEP performs 
symbolic calculation of the squared matrix element and numerical integration over multiparticle phase 
space for the physical process. The goal of the CompHEP/SUSY package is to develop the description 
of the Minimal Supersymmetric Standard Model (MSSM) and some of its extensions in the CompHEP 
format. 

We have used the description of the MSSM physical spectrum and the Lagrangian for 3 generations of 
quarks and leptons given in |2j . The LanHEP software package || was used to generate the Feynman rules 
psj . for the MSSM Lagrangian. LanHEP allows to write the Lagrangian in terms of initial fields (multiplets) 

before SU(2) x U(l) breaking, to use the superpotential formalism and the 2-component spinor notation. 
Thus the LanHEP input model description is compact and close to textbook formulae. LanHEP outputs 
O the Feynman rules as CompHEP tables describing the physical model (a LaTeX output is also possible). 

The Feynman rules were successfully compared with the ones listed in Q , as well as with many results of 
other publications. The comparision to GRACE/SUSY program was successfully done for many processes. 

We have used the notation of independent parameters same as in the ISASUSY package || which is 
more popular than one of (we give explicit mass matrices below). We also add the effective Higgs 
potential to the tree level MSSM Lagrangian to incorporate the radiative corrections to the masses and 
couplings of Higgs bosons in a gauge invariant way. 

Several CompHEP/SUSY models can be found on the WWW page 
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tittp : / / theory . sinp . msu . ru/~semenov/mssm . html 



The basic CompHEP/SUSY implementation of MSSM assumes massless fermions of the 1st and 2nd 
generation, and left /right-handed states for the corresponding sfermions. Stops, sbottoms and staus are 
mass eigenstates, the mixing due to fermion masses and trilinear soft coupling constants is taken into 
account in the calculation of the sfermion masses and interaction vertices. In the existing extension of the 
basic CompHEP/SUSY model, the fermion masses and sfermions mixing of the 2nd generations are also 
taken into account. We neglect the mixing of the different flavour left-handed squarks due to off-diagonal 
soft masses and Cabbibo-Kobayashi-Maskawa (CKM) matrix elements 0. To keep our model explicitely 
gauge invariant, we assume that only the Cabbibo angle contribute to the CKM matrix, and that the 
soft masses for the 1st and 2nd sfermion generations are equal. 

Besides the basic MSSM, there is a CompHEP/SUSY model with R-parity violation. Some other 
extensions will be described in forthcoming publications. 



2 Model parameters 

First, the user has access to the Standard Model parameters. We list those parameters with their default 
values. Some parameters, like s-, c-quark masses, CKM parameters do not appear in the basic CompHEP 
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MSSM since we neglect them. (Note that the default values might be changed in future versions). 



EE 


0.31223 


PJeotrornae'notio oonnlinj?' constant at 7V/V 

J-^JXV-' V_> IjX L X \jU LI L./XXXX &L VjUllU Lj CX I ± Li CXj \J ± VJ. £ 


GG 


1.238 


Stronc ronnlinp' constant at A//V 


SW 


0.473 


simis ot trio AA/oinnorc ancle 1 sin mi/ 


sl2 


0.221 


Parameter of OTCM matrix fsinns of (HanniHo ando^ 

_L CX X ClrXXXV-; UL'l w X v/X VlVX XXXCXj Lj X 1A \ ulll LXO VJ X V.. CX l_r l_» X CXXXtL X\_> f 


s23 


0.041 


Parameter of CKM matrix 


sl3 


0.0035 


Parameter of CKM matrix 


MZ 


91.1884 


mass of Z boson 


Mm 


0.1057 


mass of muon 


Mtau 


1.777 


mass of tau-lepton 


Mc 


1.42 


mass of e-nnarl^ 

xxxciijij vjl Vj uuciiiri 


Ms 


0.2 


mass of s-quark 


Mtop 


175 


mass of t-quark 


Mb 


4.62 


mass of b-quark 


wtop 


1.7524 


width of t-quark 


wZ 


2.4944 


Z boson width 


wW 


2.08895 


W boson width 



The following parameters are specific for the MSSM, they describe essentially MSSM extended Higgs 
sector and soft supersymmetry breaking potential. Our conventions correspond to the ISASUSY notation. 



tb 


tan/5 


Ratio of the vacuum expectation values of the Higgs doublets 


MH3 


M A 


Mass of the pseudoscalar Higgs 


mu 




Higgs mass parameter 


MG1 


Mi 


U(l) gaugino mass 


MG2 


M 2 


SU (2) gaugino mass 


MG3 


M 3 


Gluino mass 


M12 


M f 


Left-handed slepton mass, 2nd generation 


M13 




Left-handed slepton mass, 3rd generation 


Mr2 


M~„ R 


Right-handed slepton mass, 2nd generation 


Mr3 


Mf R 


Right-handed slepton mass, 3rd generation 


Mq2 


Mo 


Left-handed squark mass, 2nd generation 


Mq3 


M n 


Left-handed squark mass, 3rd generation 


Mu2 


M- CR 


Right-handed u squark mass, 2nd generation 


Mu3 


M iR 


Right-handed u squark mass, 3rd generation 


Md2 


Ms R 


Right-handed d squark mass, 2nd generation 


Md3 


M l R 


Right-handed d squark mass, 3rd generation 


Atop 


A t 


Trilinear coupling for top quark 


Ab 


A b 


Trilincar coupling for b-quark 


Atau 


A T 


Trilinear coupling for r-lepton 



We assume here that soft masses for the first generation coincide with the parameters of the second 
generation, but in future extensions the parameters Mil , Mrl, ... may appear. 

The derived parameters needed for a given process such as physical masses, mixings and so on are 
computed from those listed above and their values can be found in the CompHEP Constraints menu. 

3 Model particles 

The notation for gauge bosons are the same as in the CompHEP Standard Model: A, Z, W+, W-, G. The 
gluino is named as ~G (with mass MG3). The notation for other particles follows with more details. 
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3.1 Charginos 

Two charginos \i ■> X2 are mass eigenstates resulting from mixing of charged wino and higgsino. Com- 
pHEP names for these two particles are ~1+, ~2+. Two mixing matrices are involved to diagonalize the 
mass matrix of charginos: 

( Zmll Zm21 \ / M 2 V2M W sin/? \ / Zpll Zpl2 \ , w . 

{ Zml2 Zm22 ) { ^2M w cos/3 fx ){ Zp21 Zp22 J = ^(MCl, MC2). 

The values of the chargino masses MCI, MC2, and mixing elements Zmij, Zpij can be found in the 
CompHEP Constraints menu. 



3.2 Neutralinos 



The four neutralinos xl are named ~ol, ~o2, ~o3, ~o4. The 4x4 mixing matrix Znij is defined in the 
bino- wino- higgsino basis, so 



Z~ 1 M A fZ n = diag(W$F.l, MNE2, MNE3, MNE4), 



where 



Mm = 



Mi 
M 2 

— Mz cos (3 smOw Mz cos (3 cos 9w 

\ Mz sin (3 sin 9w — Mzsin/3cos^vK 



— Mz cos (3 sin 9w 
Mz cos (3 cos 




Mz sin /3 sin 0^ \ 
— Mzsin/3cos 







The values of Znll, ... Zn44 and masses MNE1 . . . MNE4 are in the CompHEP Constraints menu. 
Note that our definition of the mixing matrix is the inverse of that commonly used, the same is true for 
charginos and sfermions. 



3.3 Sleptons 

For each of the 3 lepton generations, there are two charged sleptons and one sneutrino. For the 3rd 
generation, the stau mixing matrix is: 

^2 _ f M l 3L +ml + M z cos2/3(-i + sin 2 6 W ) m T (A T - jUtan/3) \ 

f ~\ m T (A T - fj, tan /3) Mf B + to 2 - M§ cos 2/3 sin 2 6^ / 

The following table summarizes the notation for particles and parameters used in the basic Com- 
pHEP/SUSY model. 



Savour 


e 




T 






v T 


fcrmion 


e 


m 


1 


ne 


nm 


nl 


left scalar 


~eL 


~mL 




~ne 


~nm 


~nl 


right scalar 


~eR 


~mR 










light scalar 






~11 








heavy scalar 






~12 








left scalar mass 


MSeL 


MSmL 




MSne 


MSnmu 


MSntau 


right scalar mass 


MSeR 


MSmR 










light scalar mass 






MStaul 








heavy scalar mass 






MStau2 








M LL 






MS1LL 








M LR 






MS1LR 








M RR 






MS1RR 








Z\L 






Z133 








Zih 






Z136 








Z\R 






Z163 








Z2R 
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Z166 
MSlth 
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The parameters Mll, Mlr, and Mrr are the elements of the slepton squared mass matrix to be 
diagonalized (in GeV 2 ); the Zn,m values set the mixing of Sl,r left /right states into the si t 2 mass states: 

S L = SiZi L + S 2 Z 2 L, 
SR = SlZlR + S2Z2R, 

and 9 is the mixing angle which allows to represent the latter formulae as 

s l = s\ cos 9 — S2 sin#, 
sr = si sin 9 + S2 cos 9. 

The 9 value is given in radians. 



3.4 Squarks 

There are two squarks corresponding to each of the 6 quarks of the Standard Model. The mass matrices 
for up and down squarks of the 3rd generation are: 

/ M| £ +m t 2 +M|cos2/3(i - fsin 2 ^) m t ( A t - H cot j3) 

Mt ~ I m t (A t - fi cot (3) Mf R +m 2 + |Af| cos2/3sin 2 



/ +m 2 + Af|cos2/3(-i + isin 2 6' 1 y) m b {A b - ^tan^) \ 

M b- \ mb(A b - (i tan 0) M? + m 2 - \M% cos 2/3 sin 2 6*^ J ' 

The following table summarizes the particles and values used in the basic model. The notation is 
similar to the one used for sleptons. 



flavour 


d 


u 


s 


c 


b 


t 


fcrmion 


d 


u 


s 


c 


b 


t 


left scalar 


~dL 


~uL 


~sL 


~cL 






right scalar 


~dR 


~uR 


~sR 


~cR 






light scalar 










~bl 


~tl 


heavy scalar 










~bl 


~t2 


left scalar mass 


MSdL 


MSuL 


MSsL 


MScL 






right scalar mass 


MSdR 


MSuR 


MSsR 


MScR 






light scalar mass 










MSbotl 


MStopl 


heavy scalar mass 










MSbot2 


MStop2 


M LL 










MSbLL 


MStLL 


M LR 










MSbLR 


MStLR 


M RR 










MSbRR 


MStRR 


ZlL 










Zd33 


Zu33 


Z2L 










Zd36 


Zu36 


Z\R 










Zd63 


Zu63 


Z2R 
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Zd66 
MSbth 


Zu66 
MStth 



3.5 Higgs particles 

The Higgs particles: h°, H°, i? ± , A are named in CompHEP as h, H, H+/H-, H3 correspondingly. 
Their masses are denoted as Mh, MHH, MHc , MH3. 

It is well known that the masses of Higgs bosons in the MSSM are significantly shifted by the radiative 
corrections. There are several techniques to compute the corrections to Higgs masses. We use the effective 
potential (see || [| and references therein) to keep the explicit gauge invariance of the CompHEP model. 
The Higgs potential now is an extension of the MSSM potential: 
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V eff = m 2 ll {H 1 H+) + m 2 22 (H 2 H+) - [m 2 2 (e#itf 2 ) + h.c] 

+ \ [\(9 2 + 3' 2 ) + Ai] {H lH + f + \ [\{g 2 + g' 2 ) + A 2 ] (H 2 H+f 

+ [\(9 2 ~ S' 2 ) + As] (H 1 H+)(H 2 H+) +[-\ 9 2 + A 4 ] |eM 2 | 2 
+ (\ 5 (eH 1 H 2 ) 2 + [A 6 (fliHj) + A 7 (# 2 ff 2 *)] (eJTifl" a ) + h.c) 

The case of all A, being zero corresponds to the original tree-level MSSM potential. 

Here the parameters mf 1; m 22 can be fixed by the condition of potential minimization, and the 
parameter m\ 2 can be fixed if we choose Ma as independent variable. Then the masses of Higgs bosons 
and mixing angle a can be derived (see [|| for formulae). 

Analytical formulae for A, are given, for example, in 0. The CompHEP notation for A,; is dlhl . . . 
dlh7. 

Many authors use different approaches to calculate Higgs masses, involving diagrammatic calculation 
rather than effective potential and RGE. So, many articles and computer programs calculate only the 
CP-even Higgs masses rather than the full Ai set. This relates in particular to the FeynHiggs package. 
Note that such a prescription bounds only several of the 7 parameters, and we have to deal with A's to 
keep our model gauge invariant. 

In the case when the Higgs masses are given, we derive the values for Ai i2 .3 ; 4 assuming that Mg, Mh, 
Mh± and a are independent variables. The masses Mh, MHH, MHc and sin a sa become free parameters. 
Then we use the FeynHiggsFast package || to evaluate these parameters. This approach allows us 
to take into account the corrections to the Higgs masses and also to self-interaction vertices (see the 
discussion in ||). 

4 Width parameters 

There is a number of parameters representing particle widths. In the current version the values for 
these parameters are set by the user. These parameters can be important in case of s-channel resonant 
production of a particle. 

• wCl , wC2 — the widths of two charginos; 

• wNE2, wNE3, wNE4 — the widths of three neutralinos (except the lightest one); 

• wSG — the width of gluino; 

• wSel, wSe2, wSmul, wSmu2, wStaul, wStau2 — the widths of charged sleptons; 

• wSne , wSnmu, wSntau — the widths of sneutrinos; 

• wSul, wSu2, wSdl, wSd2, wSsl, wSs2, wScl, wSc2, wSbotl, wSbot2, wStopl, wStop2 — the 

widths of squarks; 

• wh — width of light CP-even Higgs; 

• wHh — width of heavy CP-even Higgs; 

• wHc — width of charged Higgs; 

• wH3 — width of CP-odd Higgs. 

5 SUGRA and GMSB models 

The number of independent parameters in the MSSM can be reduced in the context of the SUGRA or 
GMSB models. More specifically, the parameters of soft SUSY breaking terms (gaugino masses Mi, 
M 2 , M3, sparticle masses Ml, Mr, Mq trilinear couplings A t , A\>, A T are computed from just a few 
parameters. 
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It is possible to link CompHEP/SUSY with the ISASUSY package @, which includes SUGRA and 
GMSB models. For this, Isajet library should be installed on the computer. ISASUGRA is used to 
calculate soft SUSY breaking masses and couplings which were CompHEP independent variables. The 
masses of the particles are calculated by CompHEP from these soft parameters as before. In particular, 
masses of CP-even Higgs particles are calculated by FeynHiggsFast package, while ISASUSY gives the 
mass of the CP-odd Higgs boson. 

SUGRA option: 



CompHEP notation 


parameters 


tb 


ratio of vacuum expectation value 


suMO 


Mq - universal mass for scalars 


suMHF 


M 1 / 2 universal mass for fermions 


suAO 


Aq universal trilinear coupling 


suMU 


sign(/i) - sign of higgsino mass parameter 


suMODE 


choose the mSUGRA model 



Here suM0DE=l means default SUGRA calculation, suM0DE=2 is SUGRA with true gauge coupling 
unification at the GUT scale. 
GMSB option: 



CompHEP notation 


parameters 


tb 


ratio of vacuum expectation value 


gmLAM 


A — SUSY breaking parameter 


gmMES 


Mmes messenger mass 


gmN5 


Y5 effective number of messenger generations 


gmMU 


sign(/i) - sign of higgsino mass parameter 


gmMODE 


should be 1 (reserved for extensions) 



6 Massive fermions of the second generation 

There is an extension of the basic CompHEP/SUSY model which takes into accont the masses of the 
fermions of the second generation and the mixing of the corresponding scalar superpartners. We introduce 
the parameters for trilinear couplings Am, Ac, As for muon, c- and s-quarks. The sfermion names and 
parameters read as: 



flavour 


A* 


s 


c 


fermion 


m 


s 


c 


light scalar 


~ml 


~sl 


~cl 


heavy scalar 


~m2 


~s2 


~c2 


light scalar mass 


MSmul 


MSsl 


MScl 


heavy scalar mass 


MSmu2 


MSs2 


MSc2 


M LL 


MSmuLL 


MSsLL 


MScLL 


M LR 


MSmuLR 


MSsLR 


MScLR 


M RR 


MSmuRR 


MSsRR 


MScRR 


Z\L 


Z122 


Zd22 


Zu22 


Z2L 


Z125 


Zd25 


Zu25 


ZlR 


Z152 


Zd52 


Zu52 


Z2R 


Z155 


Zd55 


Zu55 


9 


MSmuth 


MSsth 


MScth 



G 



7 R-parity violation 



We have implemented in CompHEP the MSSM extension with R-parity violation. We incorporate here 
only the trilinear terms breaking the leptonic number. The superpotential reads: 

Wrpv — Xijk^LiLjEk + X' i: j k eLiQjDk, 

where the indices i, j, k enumerate 3 generation of matter fields. 

The CompHEP notation reads Llijk for Ay/, and Lqijk for A^- fe . A lot of processes were compared 
with the SUSYGEN package [M with good agreement. The values for Ayfc should be set twice greater 
in SUSYGEN than in CompHEP due to the different convention. 
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